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Abstract 

In a series of papers which are either published |14l I13| or available as preprints |16l 11511121117| Hadjes¬ 
fandiari and Dargush have reconsidered the linear indeterminate couple stress model. They are postulating 
a certain physically plausible split in the virtual work principle. Based on this postulate they claim that 
the second-order couple stress tensor must always be skew-symmetric. Since they use an incomplete set of 
boundary conditions in their virtual work principle their statement contains unrecoverable errors. This is 
shown by specifying their development to the isotropic case. However, their choice of constitutive parameters 
is mathematically possible and still yields a well-posed boundary value problem. 
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1 Introduction 

Among higher gradient elasticity models [551 HZl 131 [Ml UHl IM] one of the very first models considered in 
the literature is the so called indeterminate couple stress model [ilia Ea US] in which the higher gradient 
contributions only enter through gradients of the continuum rotation, i.e. the total elastic energy can be 
written as 

lT(Vw, V(Vu)) = Wiin(symVM) + Wcurv(VcurlM). 

In general, higher gradient elasticity models are used to describe mechanical structures at the micro- and 
nano-scale or to regularize certain ill-posed problems by means of these higher gradient contributions. 

In a series of papers which are either published [Mills] or available as preprints [m ESI El El] Hadjesfandiari 
and Dargush have reconsidered the linear indeterminate couple stress model. They are postulating a certain 
physically plausible split in the virtual work principle. Based on this postulate they claim that the second-order 
couple stress tensor must be skew-symmetric. Since their development has spread considerable confusion in the 
field of higher gradient elasticity, we were prompted to carefully re-examine their claim. In doing so we hope to 
contribute an important clarification in the field and to put an end to the above mentioned confusion. 

In the course of our re-examination it turned out that the boundary conditions in the classical indeterminate 
couple stress theory have never been correctly derived. In [32l [23l [7] we provide, for the first time the consistent 
and complete boundary conditions for the classical indeterminate couple stress model. In doing so, we find the 
underlying error in the argument by Hadjesfandiari and Dargush ElElEaEaEIlEZ]. While Hadjesfandiari 
and Dargush start with the linear general anisotropic couple stress response and only later specify to isotropy, 
for definiteness, we consider from the outset the linear isotropic indeterminate couple stress case. Nevertheless, 
our development is essentially independent of any isotropy assumption. It is clear that exhibiting the errors in 
their development for the simpler case of isotropy is sufficient for invalidating their claim. 

Before discussing the papers by Hadjesfandiari and Dargush [HElETlEi] we will first recall the indeter¬ 
minate couple stress model in its accepted format as far as kinematics and equilibrium equations are concerned. 
We also need to introduce the new set of traction boundary conditions which rectify the shortcomings in all 
previous papers. For comparison, the up to now accepted traction boundary conditions are also presented. 
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In the light of this new framework, we try to follow the argument given by Hadjesfandiari and Dargush 
[la na na [n] as closely as possible. We will show that their implicitly formulated requirement, recast as a 
physically plausible postulate by us, leads to the skew-symmetry of the couple stress tensor if and only if the 
classical traction boundary conditions are assumed. However, within the corrected format of traction boundary 
conditions no similar conclusion is possible. 

Despite the finally erroneous claim by Hadjesfandiari and Dargush we recognize their work in being the reason 
to reconsider the indeterminate couple stress model and to finally find the underlying error which occurred in 
the accepted boundary conditions and not directly in Hadjesfandiari and Dargush’s work. 


2 Notational agreements 

In this paper, we denote by the set of real 3x3 second order tensors, written with capital letters. 

For a,b G we let (a, &)r 3 denote the scalar product on with associated vector norm ||a||g 3 = (a, 0 )^ 3 . 
The standard Euclidean scalar product on is given by (X, F)r 3 x 3 = tr(XF’^), and thus the Frobenius 

tensor norm is ||X|p = (X, X)r 3 x 3 . In the following we omit the index The identity tensor on 

r3x3 denoted by 1, so that tr(X) = {X, 1). We adopt the usual abbreviations of Lie-algebra theory, 

i.e., so(3) := {X G = —X} is the Lie-algebra of skew symmetric tensors and s[(3) := {X G 

]^ 3 x 3 I tr(X) = 0} is the Lie-algebra of traceless tensors. For all X G we set symX = -b G Sym, 

skewX = ^{X — X"^) G so(3) and the deviatoric part devX = X — ^ tr(X)l G s[(3) and we have the orthogonal 
Cartan-decomposition of the Lie-algebra 01 ( 3 ) 

0 [( 3 ) = {sl(3) n Sym(3)} ©so(3) ©R-1, X = devsymX + skewX + - tr(X)l. (2.1) 

o 

Throughout this paper (when we do not specify else) Latin subscripts take the values 1, 2, 3. Typical conventions 
for differential operations are implied such as comma followed by a subscript to denote the partial derivative with 
respect to the corresponding cartesian coordinate. We also use the Einstein notation of the sum over repeated 
indices if not differently specified. Here, we consider the operators axl : so(3) —)> R^ and anti : R^ —)• so (3) 
through 

(axM). = A,, A„=(»lAx.. (a„ti(„))., =A«=a„ti(axlA)«, (2.2) 

for all n G R^ and A G so(3), where is the totally antisymmetric third order permutation tensor. We 
recall that for a third order tensor E and X G R^^^, u G we have the contraction operations E : X G R^, 
E. n G R^^^ and X. G R^, with the components 


(E : X)i = Eyfc Xkj , (E. v)ij = Eijk Vk , {X. v)i = Xij vj. (2.3) 

For multiplication of two matrices we will not use other specific notations, this means that for H, H G R^^^ we 
are setting = A^kBkj- 

We consider a body which occupies a bounded open set H of the three-dimensional Euclidian space R^ and 
assume that its boundary dLl is a piecewise smooth surface. An elastic material fills the domain H C R^ and 
we refer the motion of the body to rectangular axes Oxi. Let us consider an open subset F of did. Here, is a 
vector tangential to the surface F and which is orthogonal to its boundary clF, t = n x is the tangent to the 
curve cIF with respect to the orientation on F. We assume that dit is a smooth surface. Hence, there are no 
singularities of the boundary and the jump |a • := [a ■ + [a ■ v]~ = ([a]'*' — [a]“) ■ v oi a across the joining 

curve cIF arises only as consequence of possible discontinuities of the corresponding quantities which follows 
from the prescribed boundary conditions on F and dil \ F, where 


lim _ 

X e dfl \ r 
dr 


[.]+:= lim [.]. 

x G r 
or 


The usual Lebesgue spaces of square integrable functions, vector or tensor fields on il with values in R, R^ 


3 


or respectively will be denoted by L^{yi). Moreover, we introduce the standard Sobolev spaces [T1 l8l [2T| 

Hi(n) = {u e L'^{Vt) I gradw e l|w||Hi(n) := I|w|li 2 (n) + IlgradulH^^j^) , 

H(curl; n) = {v G L^{n) \ curlw € L‘^{n)}, ||w||H(curi;n) — ll^lli^cn) + l|curlr;||i 2 (o) , (2.4) 

H(div;f2) = {r; e L‘^{n)\divv G L^(f2)}, IMlH(div;n) — \Mh{n) + l|div?^||i 2 (o) , 

of functions u or vector fields v, respectively. Furthermore, we introduce their closed subspaces iJg(n), 
Ho (curl; n) as completion under the respective graph norms of the scalar valued space (7“(H) (the set of 
infinitely differentiable functions with compact support in H). For vector helds v with components in H^(H), 

i.e. V = (ui,U 2 ,U 3 )^ ,Vi G H^(H), we define Vu = ((Vui)^, (V 1 ) 2 )^, (Vwhile for tensor fields P with 
rows in H(div;H), i.e. P = [P^,P 2 ,P^), Pi G H(div;H) we define DivP = (divPi,divP 2 ,divP 3 )^. The 
corresponding Sobolev-spaces will be denoted by P^(fl), P^(Div; H). 


3 The classical indeterminate couple stress model 


We are now shortly re-deriving the classical equations based on the V[axl(skew VM)]-formulation of the inde¬ 
terminate couple stress model. This part does not contain new results, see, e.g., [53] for further details, but is 
included for setting the stage for this contribution. 

The linear isotropic indeterminate couple stress problem can be viewed as a minimization problem 



^ II symVull^ 


- [tr(Vu)]^ 




(Vcurlu) — (/, u) 


dv 


-G min. w.r.t. u, 


(3.1) 


subjected to geometric and mechanical boundary conditions, in part depending on the form of Wcurv(Vcurlu), 
which will be specified later on. 

In the following, in order to place the subject in the literature, we outline some curvature energies proposed 
in different isotropic second gradient elasticity properly models: 

• the indeterminate couple stress model (Grioli-Koiter-Mindlin-Toupin model) [9l [2l [T9l [28l l42l l40l ITO] 
in which the higher derivatives (apparently) appear only through derivatives of the infinitesimal continuum 
rotation curlu. Hence, the curvature energy has the equivalent forms 


bFcurv(Vcurl u) 


dLl 


^ II symVcurl u|p -I- ^ || skew Vcurl m|P 
[ai II sym V[axl(skew Vu)]|p -I- a-i || skew V[axl(skew Vu)] |p] 
^ II dev symVcurl u|p + ^ || skew Vcurl u|p . 


(3.2) 


We remark that the spherical part of the couple stress tensor remains indeterminate since tr( Vcurl u) = 
div(curlM) = 0. In order to prove the pointwise uniform positive definiteness it is assumed following 
m, that ai > 0 ,02 > 0. Note that pointwise uniform positivity is often assumed m when deriving 
analytical solutions for simple boundary value problems because it allows to invert the couple stress- 
curvature relation. It is clear that pointwise positive definiteness is not necessary for well-posedness m- 
Mindlin [281 p. 425] explained the relations between Toupin’s constitutive equations [H] and Grioli’s [S] 
constitutive equations and concluded that the obtained equations in the linearized theory are identical, 
since the extra constitutive parameter rj' of Grioli’s model does not explicitly appear in the equations 
of motion but enters only the boundary conditions. The same extra constitutive coefficient appears in 
Mindlin and Eshel’s version m- 


• the modified - symmetric couple stress model - the conformal model. On the other hand, in the 
conformal case [361135] one may consider that 0:2 = 0, which makes the couple stress tensor m symmetric 
and trace free. This conformal curvature case has been considered by Neff in [33], the curvature energy 
having the form 

IHcurv(Vcurl u) = fiLl ^ II symVcurl u|p = || devsymV[axl(skew Vu)|p. (3.3) 
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Indeed, there are two major reasons uncovered in [86] for using the modified couple stress model. First, in 
order to avoid singular stiffening behaviour for smaller and smaller samples in bending [34] one has to take 
0.2 = 0. Second, based on a homogenization procedure invoking an intuitively appealing natural “micro¬ 
randomness” assumption (a strong statement of microstructural isotropy) requires conformal invariance, 
which is again equivalent to 02 = 0. Such a model is still well-posed [18] leading to existence and 
uniqueness results with only one additional material length scale parameter, while it is not pointwise 
uniformly positive definite. 

• the skew-symmetric couple stress model - the non-conformal model. Hadjesfandiari and Dargush 
strongly advocate [ladiiiis] the opposite extreme case, oi = 0 and 02 > 0, i.e. they used the curvature 
energy 

IFcurv(Vcurl u) = II skew V(curlit)||^ = fj.Ll 02 || axlskew V(curlM)|p = 4 /r 02 ||curl (curlu)||^. 

In that model the nonlocal force stresses and the couple stresses are both assumed to be skew-symmetric. 
Their reasoning, based in fact on an incomplete understanding of boundary conditions is critically discussed 
in this paper and generally refuted, while mathematically it is well-posed, which will also be shown. 


3.1 Equilibrium and constitutive equations 

Taking free variations Su G in the energy kF(sym Vw, Vcurlw) = lFiin(sym Vm) -I- ITcurv(Vcurlu), where 

— ItrlV?/,!!^ = II. II devsvm V?/.ll^ -I- ^ 


Wiin(sym Vit) =fj, II sym Vu|p -1- ^ [tr(VM)]^ = M II devsym Vu|p -I- ^ ^ ^ [tr(Vw)]^, 
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Ikcurv(VcurlM) =plL\ [oi II devsym V[axl(skew Vu)]|p -|- 02 || skew V[axl(skew Vw)]!!^], 


(3.4) 


we obtain the virtual work principle 


[ WCVu + tVSu) dv 

dt Jn 


t=o 


In _ 


2^ (sym Vu, sym Vi5m)-I-A tr(VM) tr(V(5ii) 


+ fiLl[2ai (devsym V[axl(skew Vu)], devsym V[axl(skew V(5u)]) (3.5) 


2q! 2 (skew V[axl(skewVu)], skew V[axl(skewV(5u)])] -I- (/, 6u) 


dv = 0, 


where / denotes the body force density. 

The classical divergence theorem leads to 


/ (Div {a — t) + f,Su) dv — / {{a — T).n,Su) dv — / (to. n,axl(skew V<5it)) da = 0, 

Jn Jdn Jdn 


an 

the virtual power work of the surface forces 


(3.6) 


where 


CTtotai = (J-T ^ Sym(3) 

total force-stress tensor 


a = 2 fj, sym Vm -I- A tr(Vit)l G Sym(3) 

local force-stress tensor 


f = - antiDiv[TO] G so (3), 

nonlocal force-stress tensor 

(3.7) 

in = fj.L^ [ofi sym V(curl u) + 02 skew V(curl rt)] 

couple stress tensor 



= [oi dev sym V(curl u) + 02 skew V(curl m)] 

= fiLl [2ai devsym V[axl(skew Vm)] -I- 2 02 skew V[axl(skew Vm)]], 


and n is the unit outward normal vector at the surface dfl. The equilibrium equation are therefore 


Div (Ttotal + / = 0. (3.8) 

Note that the local force-stress tensor a is always symmetric, the nonlocal force-stress tensor r is automatically 
skew-symmetric, while the second order hyperstress tensor (the couple stress tensor) to may or may not be 
symmetric, depending on the material parameters. The asymmetry of force stress is a hidden constitutive 
assumption, compare to [7]. 
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3.2 Boundary conditions 

3.2.1 Classical (incomplete) Grioli-Koiter-Mindlin-Tiersten boundary conditions 

Since the variations at the boundary and the interior can be assigned independently, see (ISJl) . we must have: 


~ / (('^ — T).n, Su) da — {m. n, axl(skew V5u)) da = 0 or equivalently 

J OO J dQ 

~ / {{cr— T).n,Su) da — 2 / {rri.n,cm\Su) da = 0. 

J dCl J dfi 


(3.9) 


This suggests 6 possible independent prescriptions of mechanical boundary conditions; three for the normal 
components of the total force stress {a — T).n and three for the normal components of the couple stress tensor. 
The possible Dirichlet boundary conditions on T C seem to be the 6 conditiont0 

u = u, axl(skewVu) = u; (or equivalently curlM = 2'u;), (3.10) 

—)• at the open subset T C dfl of the boundary (3+3 boundary conditions). 


for two given functions m, w : R'^ —)• R" 
However, following Koiter we note 


is known. Then curl 


exists and 


Remark 3.1. [independent variations and curl] Assume u G and u 

for any open subset T C dkl the integral Jp(curlu, n) da is already known by Stokes theorem, while Jp(curlM, r) da 
is still free, where r is any tangential vector field on the open set T C dfl. Only the two tangential components 
of curl u may be independently prescribed on an open subset of the boundary. 

Already, Mindlin and Tiersten [5S] have rightly remarked that also in this formulation only 5 mechanical 
boundary conditions can be prescribed. They rewrote (13.91) in a further separated form 

— / {{a — T).n —-n X'V[({sym.fh).n,n)],Su) da — / ((1 — n ® n) fh.n, (1 — n 0 n) [axl(skew V(5 m)]) da 
Jon +_^ Jqq V-,- 


performs already 
work only against Su 


cannot be assigned arbitrarily independent of Su 
since it still contains certain 
’’tangential” derivatives of Su G 0°°(f^) 


(3.11) 


f {{a — T).n ——n X y[{{symm).n,n)],6u) da — 2 ( ((1 — n (g) n) m.n, (1 — n (g) n) [curldu]) da = 0, 

JdQ. 2_^ 


performs already 
work only against Su 


which already shows that the second term in (13.91) still contains contributions which perform work against 5u, 
namely x V[((sym m).n, n)], while the remaining higher order term (1 — n ® n) fh.n performs work against 
combinations of second derivatives. Mindlin and Tiersten [28] concluded that 3 boundary conditions derive from 
the first integral (correctly) and two other from the second integral, since [28l p. 432] “the normal component of 
the couple stress vector [(m.n,n) = (sym fh.n, n)] on” dO. “enters only in the combination with the force-stress 
vector shown in the coefficient of’ 5u “in the surface integral (our first term on the right hand side of (I3.1ip i. 

Therefore, Mindlin and Tiersten [55] concluded that the boundary conditions consists in: 

• Geometric boundary conditions on T C 512: 


u 

(1 — n®n). axl(skewV m) 
or (1 — n (g) n). curia 


(1 —a®n). axl(skew Va°), 
2 (1 — n (g) n). curia®, 


(3 be) 
(2 be) 


(3.12) 


for a given function a® : R® ^ R® at the boundary. The latter condition prescribes only the tangential 
component of axl(skewVa). Therefore, we may prescribe only 3+2 independent boundary conditions. 

'^as indeed proposed by Grioli in concordance with the Cosserat kinematics for independent fields of displacements and 
microrotation. 
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• Traction boundary conditions on dH. \ F: 

(ct — T).n— -n X V[((symm).n, n)] = t, 

(1 — n ® n) fh.n = (1 — n ® n) 


traction (3 be) 

“double force normal traction” (2 be) (3.13) 


for prescribed functions t, 5 : at the boundary. 

However, while (13.121) and (13.131) correctly describe the maximal number of independent boundary conditions 
in the indeterminate couple stress model and these conditions have been rederived again and again by Yang et 
al. [13], Park and Gao m, m, etc. among others they are not correct. This is explained in the following 
paragraphs. 


3.2.2 Complete (corrected) traction boundary conditions 

The indeterminate couple stress model is not simply obtained as a constraint Cosserat model [40] , i.e. assuming 
that A = axl(skewVu). In the indeterminate couple stress model the only independent kinematical degree of 
freedom is u. We understand that the indeterminate couple stress model constructed as a constraint Cosserat 
model represents at most an approximation of the indeterminate couple stress model, in the sense that the 
boundary conditions are not correctly/completely considered. We remark that the quantity {fh.n, (1 — n ® 
n)[axl(skew V(5u)]) does still contain contributions performing work against Su alone (even though there is 
a projection 1 — n 0 n involved), which can be assigned arbitrarily and are therefore somehow related to 
independent variation Su. This case is not similar to the Cosserat theory in which we assume a priori that 
displacement u and microrotation A G so(3) are independent kinematical degrees of freedom. 

At this point, it must also be considered that the tangential trace of the gradient of virtual displacement 
can be integrated by parts once again and that the surface divergence theorem can be applied to this tangential 
part of Vi5m. Let n be the unit normal vector at a considered surface point. As it is well known from differential 
geometry, the projectors TL — n^n and n®n allow to split a given vector or tensor field in one part projected on 
the plane tangent to the considered surface and one projected on the normal to such surface (see also [SlIMlIl] 
for details)). For the sake of simplicity we assume that dVL is a smooth surface of class . As in the notation 
section, we consider the curve 5r which joins the open subsets F and \ F of the boundary. Therefore, in our 
case and in our abbreviations, the surface divergence Up- 58, ex. 7] reads: 

/ Div"® (u) da := / (1 — n (g) n, V ((1 — n (g) n) • z;)]) da = / |(u, zz)|ds. (3.14) 

Jan Jan Jar 

for any field z; G R^. 

Indeed, using the surface divergence theorem, we have obtained in [33] that 


~ / {{cr — T).n,Su) da — / (ztz,. n, axl(skew Vdzt)) da 

Jan Jan 

f I 1 

= ~ / {(<7--antiDiv[TO]). n —-rz X V[((symr7z,).n,rz)] 

Jan 2 2 


(3.15) 


— -{V[(anti[(l — n (g) n)m. n]) (1 — n (g n)] : (1 — n (g n), Su) da 


- f ((1 — n (g n) anti[(l — n (g n)m. Tzj.zi, Vdzi.rz) da — - j (|anti[(l — n (g n)m. nj. zz], dzt)ds. 
2 Jan '• -V-' 2 Jgp 


completely 5u-independent second 
order normal variation of the gradient 


Hence, there are indeed two terms 

(ct — - antiDiv[m]). n — —n x V[((symTO).n, n)] — -V[(anti[(l — n (g n)fh. n])(l — n g) n)] : (1 — n g) n) 

and 

|anti[(l — rz g n)fh. n]. v\ 
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which perform work against while only the term 


(1 — n 0 n) anti[(l — n 0 n)m. n].n 

is related solely to the independent second order normal variation of the gradient \^6u.n. This split of the 
boundary condition is not the one as obtained e.g. by Gao and Park |38j and seems to be entirely new in the 
context of the indeterminate couple stress model. 

Therefore, we need to adjoin on dfl the following complete set of boundary conditions: 

• Geometric (essential) boundary conditions on T C dfl: 

u = (3 be) (3.16) 

(1 — n 0 n){\7u).n = (1 — n 0 n)(VM°).n, (2 be) 

where —>■ is a prescribed function (i.e. 3+2=5 boundary conditions), as in |28) . 

• Traction boundary conditions on dfl \ T: 

(tr — T).n — X V[((symTO).n,n)] 

— ^{V[(anti[(l — n 0 n)TO. n])(l — n 0 n)] : (1 — n 0 n) = t, (3 be) (3.17) 

(1 — n 0 n) anti[(l — n 0 n)TO. n].n = (TL — n^n)g (2 be) 

where t, g : R^ —>• are prescribed functions on 911 \ T. 

• Jump boundary conditions on ST C dfl: 

|anti[(l — n 0 n)m. n]. z/]| = tt, (3 be) (3.18) 

where tt is prescribed on ST and leads to 3 boundary conditions. 


4 The Hadjesfandiari and Dargush’s postulate 

Let us now turn to Hadjesfandiari and Dargush’s far reaching claims. In the abstract of their paper |13j 
Hadjesfandiari and Dargush write: 

“By relying on the definition of admissible boundary conditions, the principle of virtual work and 
some kinematical considerations, we establish the skew-symmetric character of the couple-stress ten¬ 
sor [m] in size-dependent continuum representations of matter. This fundamental result, which is 
independent of the material behavior [e.g. isotropy], resolves all difficulties in developing a consistent 
couple stress theory. ” 

In their appendix of [TU p. 1263-1264] they add: 

“Therefore, the present determinate theory is not mathematically a special case of [the] indeterminate 
[couple stress] theory obtained by letting [ai = 0]. This is just a coincidence for the linear isotropic 
case where equations in both theories have some similarities. It should be realized that the determinate 
theory is not simply about fixing the constitutive equations for [the] linear isotropic couple stress 
theory of elasticity. The present size-dependent couple stress theory is the consistent couple stress 
theory in continuum mechanics. This has been achieved by discovering the skew-symmetric character 
of the couple-stress tensor [m]. Mindlin, Tiersten and Koiter did not recognize the mean curvature 
tensor [skew[Vcurlu]] as the consistent measure of deformation in continuum mechanics.’^ 

■^In |14l p. 1283, A.16] Hadjesfandiari and Dargush erroneously take the strict positivity of curvature parameters ai,a 2 > 0, used 
initially by Koiter, Mindlin and others, as belonging to the definition of the indeterminate couple stress model. This is clearly not 
the case. It can be shown that both limit cases ai > 0, 02 = 0 (modified couple stress theory) and ai = 0, 02 > 0 (Hadjesfandiari 
and Dargush choice) are mathematically admitted with the provision of using the attendant correct boundary conditions being 
defined by the solution space depending on the values of ai, 02 , see Section [6] 



In our understanding this claim is completely unfounded. Assuming isotropic response, we treat their 
“consistent and determinate” couple stress theory as the linear and isotropic accepted indeterminate couple 
stress model with constitutive parameters oi = 0 , 02 > 0 in (1X71) . 

Turning to their most important claim regarding the skew-symmetry of the coule stress tensor fh we will ex¬ 
hibit their line of thought. Their reasoning is based on their fundamental hypothesis that the normal component 
of the couple stress traction vector {fh.n, n) should vanish on any bounding surface of an arbitrary volum^. In 
the case oi = 0 we have rh S so(3) and (fh.n,n) = 0 is satisfied automatically. 

In our notation the argument of Hadjesfandiari and Dargush is given as follows [El p. 12-13]: 

“From kinematics, since [w”" := (curlu,n)] is not an independent generalized degree of freedom, 
its apparent corresponding generalized force must be zero. Thus, for the normal component of the 
surface couple vector [fh.n], we must enforce the condition 

[{fh.n, n) =0]. 

[...] we notice that the energy equation ean be written for any arbitrary volume with arbitrary surface 
within the body. Therefore, for any point on any arbitrary surface with unit normal n, we must have 

[{fh.n, n) =0]. 

Since niUj is symmetric and arbitrary in [..], [m] must be skew-symmetric. Thus, [fhF = —fh] This 
is the fundamental property of the couple-stress tensor in polar continuum mechanics, which has not 
been recognized previously. ” (see also [El p. 2500]). 

Let us interpret this statement. We rephrase it for our purpose through the formulation of an implicit 

Postulate 4.1. [Hadjesfandiari and Dargush as we understand it] In any extended continuum model 
only the total force stress traction vector {a — f).n should perform work against the independent 
virtual displacement Su at the part dfl \ T 0 / the boundary dfl, where traction boundary conditions 
are applied. 

Remark 4.2. Incidentally, this postulate is automatically satisfied in classical elasticity, Cosserat and micro- 
morphic models mf . since the possible variations of the field variables are independent anyway. Whether such 
a postulate can be satisfied in a higher gradient continuum is the concern of Hadjesfandiari and Dargush. We 
will see that this is not possible. 

Hadjesfandiari and Dargush apply this postulate to the classical (incomplete) Mindlin and Tiersten’s format 
of the boundary conditions, namely (13.1211 and (I3.I3L Inspection of the indeterminate couple stress model 
within the framework of these erroneous classical boundary conditions, see e.g. (Ol) and (lOD . shows that 
choosing symm = 0 is indeed sufficient for this postulate to be satisfied. 

Hadjesfandiari and Dargush accept 


(T.n = (T„, m.n = to„, (4.1) 

and they also realize correctly that the number of geometric or mechanic boundary conditions is 5 since the 
tangential component of the test function 6u cannot independently be varied, which is by now well established. 
To this aim, they spilt the term 


IdV 


(m.n, axlskewVdu) da = 2 / {m.n, curl 6u) da (4.2) 

Jdv 

= 2 / {{n (Si n) fh.n, curls u) da-f 2 / {{1 — n®n) fh.n, curl Su) da 

JdW -V-" Jav' ---' 

normal part tangential part 

= 2 / {fh.n, {n n). curl Su) da-f 2 / {fh.n,{l — n^n). curl Su) da 

Jav '-V-" Jav ' ---' 


normal part 


tangential part 


^They also wrote |17l p.l3]: “...the corresponding generalized force must be zero and, for the normal component of the surface 
moment-traction vector [m.n], we must enforce the condition [(m.n, n) = 0].”. 
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on any arbitrary subdomain C into its tangential and normal part. They observe that the normal part 


((n ® n) m.n, curl Su) da 


IdV 



{(sym m).n, n) 

generalized force 


(curl (5u, n) 

^-V-^ 

no independent 
degree of freedom (^) 


da 


(4.3) 


cannot be prescribed independently of 5u, since indeed their w”" := (curl(5u,n) cannot independently be 
prescribed due to Stokes theorem. To avoid a somehow felt inconsistency, they state that the corresponding 
generalized force must be zero and, for the normal component of the surface moment-traction vector, they 
enforce accordingly the (misguided) condition 


((sym m).n,n) = 0 (4-4) 

generalized force 

on any arbitrary subdomain V <Z ^ having the boundary dV. 

The equilibrium equations considered by Hadjesfandiari and Dargush (ai = 01 [13] read therefore 


Div tTtotal + / = 0, 


(4.5) 


where the total force stress is given by 

Irtotai = cr — r 0 Sym(3) total force-stress tensor 

a = 2 fi sym Vu -I- A tr(Vu)l S Sym(3) local force-stress tensor 

f = i antiDiv[TO] G so(3), nonlocal force-stress tensor (4.6) 

TO = fiLla 2 skewV(curlM) = 2fiLla2 skew V[axl(skew Vu)]] G so(3) couple stress tensor. 

To the equilibrium equation, Hadjesfandiari-Dargush adjoin on the following boundary conditions 

• Geometric boundary conditions on T C dfl: 


f (1 — n 0 n). axl(skewVu) = (1 — n ® n). axl(skew Vu°), , . 

1 or (1 — n ® n). curlu = 2 (1 — n ® n). curlu*^, ^ ' 

for a given function ^ at the boundary. The latter condition prescribes only the tangen¬ 

tial component of axl(skewVu). Therefore, they prescribe(correctly) only 3-1-2 independent boundary 
conditions. 

We may consider the equivalent geometric boundary conditions 


u = uq on r, 

(1 — n (g) n) Vu. n = (1 — n (g) n) Vuq. n on T, 


(4.8) 


where mq : is given, i.e. 3-1-2 boundary conditions. 

• Traction boundary conditions on did \ T: 


The Hadjesfandiari-Dargush-choice: 


(a + T).n 
(1 — n 0 n)fh.n 


t on dil \ r, 

{1—n®n).h on 914 \ T, 


(4.9) 


where t, ft, : are given. 
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5 No reason for the Hadjesfandiari-Dargush formulation with skew- 
symmetric couple stress tensor 


Let us try to follow the argument of Hadjesfandiari and Dargush (bona fide). It follows that the split considered 
by Hadjesfandiari-Dargush is not complete, since they nowhere do use the surface divergence theorem and they 
do not prescribe fully independent geometrically boundary conditions. We see, on the contrary, that, when 
integrated over 9D, even upon the Hadjesfandiari and Dargush restriction the remaining tangential part from 
(US), namely 



n 0 n)m.n, curl 6u) da = 

tangential part 


'an 


I (V[(anti(TO. n)) (1 — n ® n)] : (1 — n (D n), Su) 
-|- ((1 — n® n) anti(m. n).n, VSu.n) > da, 

'• -V- 

completely ^^-independent second order 
normal variation of gradient 


(5.1) 


still contains some parts with independent degrees of freedom, performing work against the normal derivative 
VSu.n. 

Further on, we explain this in more detail. Looking at the Mindlin and Tiersten approach, see also [23) . in 
the framework of Hadjesfandiari-Dargush’s assumption m S so(3), we get 


(m. n, axl(skew V5u)) = ((1 — n 0 n) m.n, axl(skew V5u)) 


-|--(n,curl[ ((symm).n,n) (Jm]) — (n x V[((symm).n,n)] ,6u), 

2 '---^ 2 V---/ 


(5.2) 


zero generalized force 
according to 

Hadjesfandiari-Dargush’s assumption 


performs work against <5i4 
and should therefore vanish 
corresponding to their Postulate 


= ((1 — n 0 n) m.n, axl(skew V(5 m)) = — n 0 n) m.n, curl^u). 


(5.3) 


it is 


Therefore ([3.lip can immediately be written as 

— / {<{(T — T).n,5u)da— / ((1 — n 0 n) m.n, (1 — n 0 n) [axl(skew Vi5m)]) da = 0, 

JdQ. '■ Jdn 

which is used and accepted by Hadjesfandiari and Dargush, see ( 1121 ). We comprehend their curvature parameter 
choice (based on the incomplete format of the independent boundary conditions): normal tractions would be 
automatically completely separated into pure total force-stress tractions and pure couple stress tractions. 

However, after integration, the second term from (15.31) ((1 — n 0 n) m.n, (1 — n 0 n) [axl(skew Vdu)]) leads 
to two new quantities: one performing work against the normal derivative S/Su.n and one still performing work 
against 5u. This means that the assumption of Hadjesfandiari and Dargush does not remove all quantities 
which may also perform work against 5u, besides the total force stress tensor (cr — T).n. This fact follows as 
a direct consequence of (|3.15l) . since, also in the Hadjesfandiari-Dargush formulation similar to 
possible for us to use the surface divergence theorem and to obtain 

~ / {{cr— T).n,6u) da — / (m. n, axl(skew Vi5m)) da 

J dQ J dfl 

1 ((cr — T).n — - V[(anti[(l — n 0 n)m. n])(l — n 0 n)] : (1 — n 0 n), Su) da 

JdQ 2 

((1 — n 0 n) anti[(l — n 0 n)fh. n].n, VSu.n) da 

i 

(|anti((l — n 0 n)m. n)]|. v, Su) ds, 

[ ((cr — r). n —-V[(anti[m. n])(l — n 0 n)] : (1 — n 0 n), (5 m) da 
Jdn 2 

- I ((1 — n 0 n) anti[m. n].n, Vdit.n) da — - j (|anti(m. n)|. z/, du) ds, 

2 Jdn 2 Jar 


IdQ 

1 

2 Jan 
1 

2 Jar 


(5.4) 
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for all variations 5u G where we have used that (m. n,n) =0 implies 

(1 — n ® n)fh. n = rh.n — n {fh. n, n) = fh. n. 


We consider this last representation as the only correct form following the Hadjesfandiari-Dargush assumption 
(to. n,n) = 0 made in (lO)) . 

Similar to the correction to the Mindlin and Tiersten’s approach with the specification that fh G so(3) (which 
is now seen to be not necessary but possible), we arrive rather, upon the Hadjesfandiari-Dargush assumption 
(to. n, n) = 0, at the traction boundary condition 


[(cr — T).n— iV[(anti[(l — n® n)m. n])](l — n®n)] : {1 — n®n) (x) 

[(1 — n®n) anti[(l — n ® n)fh. n\.n\ (x) 
on on \ r, while on 9r we have to prescribe the jump conditions 

{|anti[(l — n® n)m. n]]|. i'}{x) 


g{x), (3 be) 

[(1 — n ® n) s] (x) (2 be) 

7r{x), (3 be) 


From (j5.5L we see finally that their intended response is not satisfied: 


(5.5) 


Remark 5.1. Assuming that {{symm).n,n) = 0 (the Hadjesfandiari-Dargush assumption) implies only that 


ion 


(V[(anti[(n ® n)m. n])(l — n®n)] : (1 


However, from 


V[(anti[TO.n])(l —n®n)] : (1 


there remains another contribution 


n ® n), Su) da = 0. 
n ® n), 


(5.6) 


V[(anti[(l — n® n)m. u])(l — n®n)] : {1 — n®n) 

which still performs work against 6u. The raison d’etre of the Hadjesfandiari-Dargush formulation, as we 
understand it, was to avoid that any parts related to “couple stress normal traction” fh.n, other than the 
contributions to the total force-stresses, i.e. other than a — t, would perform work against Su. Our calculation 
shows that this requirement can never be satisfied in any higher gradient elasticity theory. 

Continuing, Hadjesfandiari tries to support the claim of a skew-symmetric couple stress tensor with some 
independent motivations in m- There, he defines the “torsion-tensor” X and the “mean curvature tensor” w 


X(u) = symVcurlu, uj{u) = skewVcurlu, (5-7) 

respectively. In Section 3.2. of [ISJ eq. 32] Hadjesfandiari claims that an inhomogeneous state u of constant 
torsional deformation X(u) = 0 cannot exist. However, the conformal mapping, see Appendix 1X1 

(j)c{x) = ^ ^2(axllF, x) X — axlIF)||x|p^ -\-\pl. A\.x -\-b, (5.8) 

where W,A G so(3), 6 G p G R are arbitrary but constant is on the one hand inhomogeneous in the 
displacement u = (jc — x but on the other hand gives precisely 


X(u) = X.{(j)c{x) — x) = sym Vcurl(/)c(x) = 0. 


(5.9) 


For instance, a simple conformal displacement field is given by 
4>c{x) = 


(x?-bxi-bxi) \ , 

( xl- (xl+xl) \ . 

( 2 X 1 

— 2X2 

-2X3 

2xiX2 = 

2xiX2 ^ V(/>c(x) = 

2X2 

2 Xl 

0 

2xiX3 / ' 

\ 2 X 1 X 3 I 

(2X3 

0 

2 Xl 


/ 0 

— 2 X2 

-2X3 \ 

/ 

0 \ 

skewV(()c(x) = 2x2 

0 

0 

^ axlskewV(()c(x) = 

-2X3 

\ 2X3 

0 

0 / 

\ 

V 2X2 / 


0 

cnr\ 4>c{x) = ( -X3 

X2 


Vcurl(()c(x) = 



sym Vcurl(/)c(x) = 0. 


(5.10) 
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6 Existence and uniqueness of the solution in the Hadjesfandiari- 
Dargush formulation 


In view of the previous discussion, we do not think that the Hadjesfandiari-Dargush formulation has a sound 
physical motivation. Nevertheless, mathematically it is possible to consider the parameter choice inherent in 
the Hadjesfandiari-Dargush formulation. Let us consider for simplicity null boundary conditions. Hence, in the 
following we study the existence of the solution in the space 

Xo = {u€ HQ{n)\cuilu € H{cuil;U), (1 — n ® n). curlu|p = O}. (6.1) 

On Xq we define the norm 

ll«b„ = (l|V«|li.(a) -I- IlskewVaxl(skewVM)||^2(Q) )' = (^l|Vu||^2(n) + J||curlcurlu||2,(j.j)^ , (6.2) 

and the bilinear form 


iin,v))= [ [: 

Jn . 


2/r (sym Vu, sym Vu) -I- A tr(Vu) tr(Vu) 

+ 2fiLla3 (skew[V axl(skew Vu)], skew[V axl(skew Vu)]) 


dv 


(6.3) 


/n . 


2/r (sym Vu, sym Vu) -I- Atr(Vu) tr(Vu) -I- 2 ^L^q ;3 (curl curl u, curl curl w) 


dv, 


where u,v € Xq. Let us define the linear operator Z ^ Rj describing the influence of external loads, 

l{v) = j^{f, v) dv for all w € Xq. We say that w is a weak solution of the problem ("P) if and only if 


{{u,v)) = l{v) for all v € Xq. 


(6.4) 


A classical solution u G Xq of the problem (P) is also a weak solution. 
Theorem 6.1. Assume that 


i) the constitutive coefficients satisfy ^ > 0, 3 A -I- 2/r > 0, 03 > 0; 

a) the loads satisfy the regularity condition f G L^(D). 

Then there exists one and only one solution of the problem (lOl) . 

Proof. In the case 03 = 0 we have the boundary value problem from classical elasticity. Further, we consider 
the case as > 0 . The Cauchy-Schwarz inequality, the inequalities (a ± 5)^ < 2(a^ + 6 ^) and the assumption 
upon the constitutive coefficients lead to 


{{u,v)) < , 

which means that ((•, •)) is bounded. On the other hand, we have 


(6.5) 


(0 


L,u)) = [ ; 
Jn _ 


2pL II sym Vu||^ 


A [tr(VM)]^ + 2 ^ 03 II skew[V axl(skew Vu)] |p] 


dv. 


( 6 . 6 ) 


for all u G Xq. Moreover, as a consequence of the properties i) of the constitutive coefficients we have that there 
exist the positive constant c 


((u,u))>c J ^11 sym Vu|p-I-II skew[Vaxl(skew Vu)]|p^ du. 
From linearized elasticity we have Korn’s inequality [53], that is 


I|Vu||l 2 (q) < C7|| sym Vm||l 2 (o) , 


(6.7) 


( 6 . 8 ) 
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skewVii -A E so (3), degenerate Cosserat model 
(j 2 // sym Vii + 2 /ic skew(Vu — A) ^ Sym(3) 
m = 2 fi skewV axl(A) E so (3) 

£ fi ||symVu|p + ^dl skew Vu — A|p + ||skewV axl(A)|p 

invariant under: u u + W.x + b, b £ 

Ah^A + W, Ty€so(3) 
well-posed: u £ axl(A) £ H(cur\]£l), 

degenerate curvature energy, 


3 -{- 2 = 5 geometric be: 
3 -j- 2 = 5 traction be: cr. 


«|p, (axl(v4), Ta>|j 


Ir 

lan\r’ 


(m.n.Ta) 


Vti = A 
constrain 

flc ^ OO 


Hadjesfandiari-Dargush model ' 

<7 2 sym Vii £ Sym(3) 

ih = 2 fi skewV axl(skewVii) € so(3) 

T = 2fiL^ anti[Div(skewVaxl(skewVw))] € so(3) 
£ ^ ^ ||symVtx|p + ||skewVaxl(skewVti)|p 
invariant under: u\-^ u-\- W.x -h 6, 6 E R^ 

Vu !->■ Vu + VF, W E so(3) 
well-posed: u S curl u G H (curl; f!) 

5 geometric be: ulp, (curlu,T q)| p 
5 traction be > 


Figure 1: A possibility of lifting the 4th.-order indeterminate couple stress model to a 2nd.-order micromorphic 
or Cosserat-type formulation formulation. Here Ta, a = 1,2 denote two independent tangential vectors on the 
boundary. 


for all functions u £ HQ(r2;r) with some constants C > 0, for bounding the deformation of an elastic medium 
in terms of the symmetric strains. Hence, using the Korn’s inequality (|6.8I) . it results that there is a positive 
constant C such that 


{u,u)>c J ^||Vu|p-I-II skew[Vaxl(skew Vm)]|P^ dx) = c ||m|| 


(6.9) 


Hence our bilinear form ((•,•)) is coercive. The Cauchy-Schwarz inequality and the Poincare-inequality imply 
that the linear operator ((■) is bounded. By the Lax-Milgram theorem it follows that (16.4|) has one and only 
one solution and the proof is complete. □ 

Remark 6.2. The Lax-Milgram theorem used in the proof of the previous theorem also offers a continuous 
dependence result on the load f. Moreover, the weak solution u minimizes on Xq the energy functional 


Hu) = f : 

Jn . 


2 pL II sym Vu|p -I- A [tr(VM)]^ -I- 2 as || skew(Curl (sym Vit))|p — (/, u) 


dv. 


7 The constrained Cosserat formulation of the Hadjesfandiari-Dargush 
model: well posedness of a degenerate Cosserat model 

Similarly to the classical indeterminate couple stress model, also the Hadjesfandiari-Dargush formulation can 
be obtained as a constrained Cosserat model. We only need to adapt the curvature energy. We consider the 
replacement skew Vxt i—>■ A £ so(3), and we obtain the energy 

W = 2 p ||symVit||^ -I- A [tr(Vxt)]^ -|- pc || skew Vm — A|p -|- 2/iL^ ||skewV axl(A)|p (7.1) 

= 2 p ||symVu|p -I- A [tr(Vxt)]^ -I- pc || skew Vm — A\\^ -\-2 pL\ ||curl axl(A)|p 

for the Cosserat model. 

Using the usual procedure, it follows that there exists a unique solution {u, A) of the corresponding mini¬ 
mization problem, i.e. to find the minimum of the energy 


I{u) = 


2 p ||symVM|p -I- A [tr(Vii)]^ + Mc|| skew Vm — A\\'^ 2 pL^ ||curl axl(A)||^ 


- {f,u) - (axl(911),axl(A)) 


dv. 


(7.2) 


where / : D —)• and 911 : D are prescribed, such that u £ i7o(D) and axl(A) £ i7(Curl;D), 

axl(A) X n|p = 0. 
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8 Conclusion 


First, Hadjesfandiari and Dargush reject the notion of additional independent degrees of freedom, which is 
purely arbitrary m p. 2496]. The “consistent” theory proposed by Hadjesfandiari and Dargush “resolving all 
difficulties...”, taking only the skew-symmetric part of the curvature tensor Vcurl u, is simply a special case 
corresponding to some vanishing moduli in the theory. As we have seen, this choice is a restriction, but not a 
necessity. 

In summary, Hadjesfandiari and Dargush have raised [181 p. 17] three major concerns in the indeterminate 
couple stress model (see also my- 

1) The body-couple is present in the constitutive relations for the [total] force-stress tensor in the 
[indeterminate couple stress] theory. 

2) The spherical part of the couple-stress tensor is indeterminate, because the curvature tensor 
[k = Vcurl u] is deviatoric. 

3) The boundary conditions are inconsistent, because the normal component of moment traction 
[(n,m.n)] appears in the formulation [of the force stress tensor]. 

Regarding the above three “serious inconsistencies” presented in m p. 17] we may answer: 

1) This is of course not inconsistent. It is well-known that the Cauchy-like total force stress tensor 
cr — r is not the constitutive stress. The constitutive stresses and couples are those arising 
from the virtual work principle (13.51) by fixing sym Vw and Vcurl rt as independent constitutive 
quantitiefl More precisely, the constitutively dependent quantities are the energetic conjugates 
of ^symVw and ^Vcurlrt, respectively. We have to note that the constitutive dependent 
quantities are not the energetic conjugates of -^u and ^curlu, respectively. 

The relation between the Cauchy-like stress and constitutive stress indeed involve the volume 
simple double and triple forces (see [B] in French, first part on the classical theory). Note that 
in the Hadjesfandiari and Dargush-formulation (14.91) . the boundary conditions would acquire 
the same form and meaning as in the classical format: the total force stress tensor cr — r would 
be the Cauchy stress tensor and the curvature of the surface normal would not intervene in the 
traction boundary conditions. However, the incorrect boundary conditions used by them (14.91) 
are not any more in the completely independent decomposition form (15.51) of the boundary 
conditions; 

2) The indeterminacy of the spherical part of the couple stress is not inconsistent. Like the pressure 
in an incompressible body, it is indeterminate in the local constitutive law but can be found 
from the boundary conditions after solving the equilibrium equations. It is a reaction stress, 
as it is well-known in the theory of continua with internal constraints. Should we say that the 
theory of incompressible bodies is “inconsistent”? Surely not, it is mathematically clear even 
though it may induce computational difficulties. 

3) As we have shown in this paper, see also [53], the boundary conditions used by Hadjesfandiari 
and Dargush m are incomplete, therefore no inconsistency occurs, see again (1531) VS. dHI). 

Finally, we like to mention that the extra constitutive parameter 77 ' of Grioli’s model does not intervene 
in the field partial differential equations. It will arise through the boundary conditions. We do not see which 
mechanical principle is violated by that. In summary, there is a fully consistent version of the indeterminate- 
couple stress model with only 3 constitutive parameters. This is the modified couple stress model, see [Ml [7] 
and a novel variant of it is recently discussed in [7] . 
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A Conformal invariance of the curvature energy and group theoretic 
arguments in favor of the modified couple stress theory 

This section is taken from [T| and included here for this contribution to be rather self-contained. An infinitesimal conformal mapping 
|33l I36| preserves (to first order) angles and shapes of infinitesimal figures. The included inhomogeneity is therefore only a global 
feature of the mapping. There is locally no shear-type deformation. Therefore it seems natural to require that the second gradient 
model should not ascribe energy to such deformation modes. 

A map 0c : —>• 1^^ is infinitesimal conformal if and only if its Jacobian satisfies pointwise V0c(a^) € M ■ 1 -|-so(3), where 

R • 1 -|-so(3) is the conformal Lie-algebra. This implies [3311361 (3^ the representation 


(pdx) = — (2(axl W, x) X — axl W||x|p) A [pi A A].x A b, 

where W, A € sci(3), b € R^, p C R are arbitrary given constants. For the infinitesimal conformal mapping 0c we note 


V0c(aj) = [(axl W, x) -h 3 1 + anti(W. x) -j- A, 
skewV0c(3^) = anti(W.3^) A, 
dev sym V0c(a^) = 0, 
sym Vcurl 0c(a:) =0, 


div0c(x) = tr[V0c(a^)] = 3 [(axl VF,^) -h^, 
symV0c(a;) = [(axlVF,^;) I? 
Vcurl0c(x) = 2W € so(3), 
skew Vcurl 0c(a:) = 2W. 


(A.l) 


(A.2) 


These relations are easily established. By conformal invariance of the curvature energy term we mean that the curvature energy 
vanishes on infinitesimal conformal mappings. This is equivalent to 

Wcurv{L)^(pc) = 0 for all conformal maps 0c, (A.3) 

or in terms of the second order couple stress tensor m := Dvcurltf^curv(Vcurlu.), 

m(Z)^0c) = 0 for all conformal maps 0c. (A.4) 

The classical linear elastic energy still ascribes energy to such a deformation mode, but only related to the bulk modulus, i.e., 

2 2 A 3 X 


^lin(V0c) = p II devsym V0c||^ + ^ ' [ti'(V0c)]^ = --- [tr(V0c)]^. 


(A.5) 


In case of a classical infinitesimal perfect plasticity formulation with von Mises deviatoric flow rule, conformal mappings are precisely 
those inhomogeneous mappings, that do not lead to plastic flow m, since the deviatoric stresses remain zero: devsym V0c = 0. 
In that perspective 


conformal mappings are ideally elastic transformations and should not lead to moment stresses. 


The underlying additional invariance property of the modified couple stress theory is precisely conformal invariance. In the 
modified couple stress model, these deformations are free of size-effects, while e.g. the Hadjesfandiari and Dargush choice would 
describe size-effects. In other words, the generated couple stress tensor m in the modified couple stress model is zero for this 
inhomogeneous deformation mode, while in the Hadjesfandiari and Dargush choice in is constant and skew-symmetri(0. 


^This observation is a further development in understanding why the Hadjesfandiari and Dargush [1211131 IT7| choice is rather 
meaningless, while mathematically not forbidden. 
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